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ABSTRACT 

Large-scale inhomogeneities and anisotropics are modeled using the Long Wavelength 
Iteration Scheme. In this scheme solutions are obtained as expansions in spatial gradients, 



Q \ which are taken to be small. It is shown that the choice of foliation for spacetime can make 
■ 

O ' the iteration scheme more effective in two respects: (i) the shift vector can be chosen so 

in : 

O ■ as to dilute the effect of anisotropy on the late-time value of the extrinsic curvature of the 



spacelike hypersurfaces of the foliation; and (ii) pure gauge solutions present in a similar 
calculation using the synchronous gauge vanish when the spacelike hypersurfaces have 
extrinsic curvature with constant trace. We furthermore verify the main conclusion of the 



^ ! synchronous gauge calculation which is large-scale inhomogeneity decays if the matter — 
c3 \ considered to be that of a perfect-fiuid with a barotropic equation of state — violates the 
strong-energy condition. Finally, we obtain the solution for the lapse function and discuss 
its late-time behaviour. It is found that the lapse function is well-behaved when the matter 
violates the strong energy condition. 

PACS numbers: 04.50.+h, 98.80.H 
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1. Introduction 

The long wavelength iteration scheme [1,2,3] has been used to model analytically 
cosmologies with large scale inhomogeneity. This scheme was first introduced by Lifschitz 
and Khalatnikov (see [4] and references therein). A similar scheme has been used by 
Salopek, Stewart, and collaborators [5] in their analysis of the Hamilton- Jacobi equation 
for general relativity. It has also been applied to Brans-Dicke theory [6] and the more 
general scalar-tensor gravity [7] . The essential part of the long wavelength iteration scheme 
is spatial gradients are considered to be small. Therefore, the scheme immediately suff'ers 
from the fact that there is no absolute space, and hence no absolute notion of spatial 
gradient, large or small. We circumvent this difficulty by using the 3-1-1 formalism to 
perform our calculations (see York [8] for an excellent review). The technical problems 
associated with the lack of an absolute space are overcome by introducing a foliation of 
spacetime, and hence a notion of time and the fiow of time. But furthermore, we will see 
that the 3-1-1 formalism can also make the long wavelength iteration scheme more effective 
at producing inhomogeneous and anisotropic cosmological models. 

Refs. [1,2,3] solve a first-order (i.e., one time derivative) form of the field equations 
where the synchronous gauge is imposed. In the language of the 3+1 formalism, the 
synchronous gauge sets the shift vector to zero and the lapse function to one. But the 
immediate simplification this gauge choice obtains, in making the field equations less com- 
plicated by having fewer fields to worry about, does not necessarily make the calculations 
less complicated (see, for instance, Piran [9], Smarr and York [10], or Hawking and Ellis 
[11]). For instance, it will be shown below that the shift vector can be used to exponentially 
damp (in time) the influence of anisotropy on the value of the extrinsic curvature of the 
spacelike leaves of the foliation. This can be significant because in the synchronous gauge 
approach obtaining the three-metric of the spacelike slices of the foliation when anisotropy 
is present is complicated and was only treated as a linear perturbation during the iterative 
part of the calculations. 

The most important difference between the formulation used here and that of Ref. 
[1] is in the choice of lapse function: we will use a lapse function such that the trace of 
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the extrinsic curvature, K, of each spacehke shce associated with the fohation is constant. 
The value of K will change from one slice to the next, and is thus a function of time, 
i.e., K = k{t). (The trace can even be taken to be the time itself [12], which is known as 
'York time.') The K = k{t) slicing condition is used because it ensures better than the 
synchronous gauge that calculated inhomogeneity is physical [12,13,14], and not merely an 
artifact of embedding unnecessarily distorted spacelike slices into an otherwise homoge- 
neous spacetime. Indeed, it will be seen below that pure gauge solutions appearing in the 
results of Comer et al. [1] do not show up here. However, we do find their other solutions 
and thus confirm the physical origin of the inhomogeneity they calculate. Furthermore, 
we also confirm their main conclusion, which is long wavelength inhomogeneity grows or 
decays, depending on the choice of the equation of state for the matter. In particular, the 
inhomogeneity decays when the matter is 'inflationary' (i.e., it violates the strong energy 
condition [11]). 

It must be stated now that we do not solve the 3+1 form of the field equations as an 
initial-value problem. That is, no explicit attempt is made to solve first the constraints — 
so as to isolate the true degrees of freedom, whose initial values make up the initial data 
set — and then use the remaining field equations to evolve the initial data set forward in 
time. Here expansions in terms of spatial gradients of a "seed" metric are made for the 
three-metric of the spacelike slices and matter variables. The zeroth-order term in an 
expansion contains no spatial gradients, the second-order term, which is next, contains 
two spatial derivatives, the fourth-order has four, and so on. (We shall demonstrate that 
the lapse function, but not the shift vector, must also be expanded in this way.) The 
coefficients in the expansions are constructed so that the constraints, as well as the other 
equations, are satisfied order- by-order. In this context, the initial- value problem becomes 
just a different way of solving the field equations; however, we will discuss it briefly in the 
concluding section. 

Building a slicing of spacetime based on keeping the trace of the extrinsic curvature 
constant has now a long history, beginning with the maximal slicing condition of Lich- 
nerowicz [15]. Since then much work has been done to determine when the K = k{t) 
condition may be used, for instance, in the context of homogeneous cosmology (see Ryan 
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and Shepley [16] for a review) or vacuum, asymptotically flat spacetimes with no restric- 
tion to homogeneity [17]. Goddard [18] has even shown for closed universes that if one 
K = k{t) slice exists, and the strong energy condition is satisfied, then the K = k{t) 
condition may be applied through the whole spacetime. More recently, Goldwirth and 
Piran [19] demonstrated numerically that the K = k{t) slicing condition breaks down for 
closed, inhomogeneous universes having 'inflationary' matter. What they find is one can- 
not prescribe a k{t) which is a monotonic function of time and simultaneously construct 
an acceptable lapse function. However, they were able to use K = k{t) for open universes, 
even during infiation, with no problems. As for the present work, there are no difficulties 
with using the K = k(t) slicing condition when the matter is 'inflationary.' We can specify 
a monotonic k{t) and still get an acceptable solution for the lapse function (in the sense 
that it remains positive for all time). On the other hand, when the matter satisfies the 
strong energy condition, caution must be exercised in the use oi K = k{t) ^ especially when 
the higher-order terms in the lapse function become comparable to the zeroth-order term. 

In Sec. 2 we give the 3-|-l form of the field equations and set the bulk of the notation. 
In Sec. 3 we discuss in general terms (ie., for arbitrary lapse function and shift vector) 
the long wavelength approximation [20]. In Sec. 4, we discuss the zeroth-, second-, and 
fourth-order equations and solutions for K = k{t) hypersurfaces. In Sec. 5, we write the 
solution for the lapse function, good to fourth-order in spatial gradients, and discuss how it 
behaves at late times. In Sec. 6, we offer some concluding remarks. Finally, some formulas 
used to construct the second- and fourth-order equations are listed in Appendix A and the 
method through which the zeroth-order solution for the extrinsic curvature is obtained is 
given in Appendix B. We will use units throughout such that G = c = 1 as well as 'MTW 
conventions (see Ref. [21]). 

2. The 3-|-l Decomposition and Field Equations 

The notion of 'time' in the 3-1-1 formalism enters via a scalar function t = t(x^), n = 
0, 1, 2, 3, whose level surfaces are spacelike. The fiow of time is introduced via a vector t^, 
which can be either timelike or spacelike. The vector which is perpendicular to the level 
surfaces of t{x^) is = —NV/j.t, where the lapse function N ensures that n'^ni, = — 1. 
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The vector is not parallel to but has the form 

r = (±;; -n^n^) t^' = N'' + Nn"" , (1) 

where -LJ1= 5^ +71^71^ is the 'projection-operator' that projects spacetime tensors into the 
level surfaces obtained from t{x^) and is the shift vector [ui^N^ — 0). From Eq. (1), it 
is seen that N must always remain positive, otherwise there would exist spacetime points 
where ti^ would be parallel to the spacelike leaves of the foliation [10]. It must also be 
emphasized that the lapse function and shift vector are arbitrary in the sense that the four 
independent functions associated with them are not determined by the field equations but 
instead must be specified by hand. 

From here on out it will be assumed that a local coordinate system = (a;°,a;*) 
exists, where the a;*, z = 1,2,3, are the coordinates of the spacelike slices associated with 
the foliation and we take t = x°. With respect to these coordinates A'"'^ = (0, N^), and the 
line interval can be written 

ds^ = -{N^- N'Ni)dt^ + 2Nidx'dt + ^ijdx'dx^ , (2) 

where = gij and Ni — jijN^ . The spatial covariant derivative compatible with 7^- is 
denoted 

The 3+1 field equations consist of four constraint equations and twelve evolution 
equations (see Ref. [8] for more discussion). They contain terms with only one time 
derivative because jij (where a dot overscript '"' means d/dt) is replaced by the extrinsic 
curvature, whose spatial components Kij are obtained from [22] 



K = — - ±'^±'^ V T? = - 

J^ij — 2 « J Vcr'tr 2 



DiNj + DjNi-% . (3) 



The inverted form of Eq. (3) thus represents six of the evolution equations (for 7^). The 
other evolution equations, and the constraints, result from projecting the free indices of 
the Einstein equations into and out of the spacelike slices of the foliation. The three inde- 
pendent projections of the stress-energy tensor T^;^ define an energy density p = n'^n^T^^, 
a three-current jTj = — n'' T^aj and spatial stress tensor Sij =1.'^ -Lj T^r- The 3-1-1 
equations are 

167rp = R- K^K{ + , (4) 
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SttJ, = Dj [kI - K51) , (5) 
SttSI = Ri + N-"- [K^DiN^ - K^DkN^^ + KKl - N'^K'^- 

N-^bib^N + N-^N^bkKi - 47r(p - S)5l , (6) 

and 

jij = biNj + bjNi - 2NKij , (7) 

where Rij is the three-dimensional Ricci tensor formed from 7^^, R = ^''Rij is its associated 
Ricci scalar, and S = Y^Sij. Note that Eq. (6) is the trace free form of -Lf-LJ G^rr = ^T^Sij, 
the removed trace being 

k = 4ttN {p + S)+ N'biK + NKlK] - D'DiN . (8) 

Lastly, the evolution equation for the determinant of the three-metric 7 (which will be 
used below) is 

7 = 27 (^AAT* - ATk) . (9) 

The field equations for a perfect fluid are VjyT^ = 0. The independent projections of 
the free index gives 

p + NbiJ" = N (k]SI + Kp) - 2Tb,N + N'bip (10) 

and 

T + NbjS^^ = N {2Kjj^ + KJ') - S'WjN - pD'N + NWjJ' - JWjN' , (11) 

where p, J'^ and Sij are obtained from the projections of the perfect fluid stress-energy 
tensor 

V = (P* + V*) u^u^ + v*g^,y . (12) 

The fluid four- velocity u^^ has the normalization u^u^ = — 1 and the pressure p* is related 
to the energy density p* via the equation of state p* = (F — 1) p*, where < F < 2. 
Letting represent the projection of the four- velocity into the slices of the foliation, i.e., 
■u* =-L^ u"', and |up = ttj, then 

p = p*(l + F|up) , (13) 
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Ji = ±rp*y^T+W , (14) 

and 

Sij = Tp* Ui uj + (r - 1) p*7y . (15) 

3. The Long Wavelength Approximation 

In this section we invoke the long wavelength approximation, where spatial gradients 
are neglected entirely. These solutions are accurate representations of cosmologies wherein 
deviation from perfect homogeneity occurs on scales bigger than the co-moving Hubble 
length (see Tomita [20], and also Ref. [1], for more discussion). It will be seen that one 
can extract a significant amount of information about the solutions before any assump- 
tions about the lapse function are made and without setting the shift vector to zero. In 
particular, since the lapse function and shift vector are specified freely, there is no a priori 
reason for their spatial gradients to vanish. 

The long wavelength equations are 

IGttpo ~ ^^Ki + , (16) 

Stt f« , (17) 

and 

N-^DiD^N - 47r(po - ^°^5)5^' , (18) 

where a '0' left-superscript (or right-subscript) on any quantity signifies it is zeroth-order in 
spatial gradients. The equations for ■jij and 7 do not change their form from the relations 
given in the last section; however, that for K does change to 

^ 47tN (po + + N ^°^Ki - DWiN . (19) 

The stress-energy divergence equations become, by forcing Eq. (17) to hold on each space- 
like slice (which means taking ^^^Ui at all times), 

Po^N ('^Si + po) (20) 
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and 

Of« ((°)5|+po5^)^iiV. (21) 

Two immediate conclusions follow from Eqs. (20) and (21): using the evolution 
equation for 7, then it can be obtained from Eq. (20) that 



po ^°^7^/^ C exp r / di DiN' , (22) 



where C depends only on x*. Eq. (21) says that DiN ^ whenever the combination 
(0)5^ + pq51 is not zero. (For a perfect fluid with the equation of state used here, the 
combination is zero only when T = 0.) Hence, even though the lapse function is in principle 
freely specifiable, consistency of the long wavelength approximation demands that it have 
a small spatial gradient. 

The remaining equations (18) and (7) are to be solved, respectively, to obtain ^'^^Kj 
and then ^'^^jij- Finding ^^^K^! requires a rewriting of Eq. (18). This is accomplished 
by using Eq. (16), and the fact that DiN and the three- velocity are both zero at this 
order. Also letting (o)/C| = U^^K^ - | ^^^KSf] , which is the trace-free part of ^^^Kj 

(modulo the •\/(°)7 factor), then Eq. (18) becomes 

(°)a:| -£n^°^/C^' ~0 , (23) 

where C^^^^ICl is the Lie derivative with respect to N^, i.e., ignoring spatial gradients of 
the extrinsic curvature and the three-metric, 

£n^°^/C^' ~ {^DkN''^ (o)/C^' + (sjDiN'' - SfDiN^^ (°)/Ci . (24) 

In order to solve Eq. (23) , it is convenient to interpret it as a linear vector equation (with 
the vector K being formed out of the components of ^^^JCj), i.e., 

K - (^DkN''^ K - MK , (25) 

where 

K = [(°)/c} (°)/c^ (o)/c? (o)/c2 (o)/c2 (o)/ci (°)/ci]^ (26) 



and 
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(27) 

with ai = £>l7V^ a2 = £'2^^^ 03 = £'3^^^ 04 = as = D2N^, ae = DsN^, 

a-j = DiN^, as = D2N^, and ag = DsN^. One can quickly verify that TrM = 0. 
Therefore, the eigenvalues of M must add up to zero. 

Even though Eq. (25) is linear in the extrinsic curvature, solving it for a general 
shift vector is a daunting linear algebra eigenvalue problem, further complicated by the 
time dependence in the matrix M. We will simplify the problem by using the freedom to 
specify, by hand, the shift vector and assume that its spatial gradient is time-independent. 
Nevertheless, finding the eigenvalues of the matrix M still means solving a ninth-order 
polynomial equation. Fortunately, a lot of progress can be made even without knowing 
explicitly the eigenvalues. By following the procedure outlined in Appendix B, then it can 
be shown that 



= exp 



(28) 



DkN^) t 

The tensor Qj {x'',t) is the analog of the so-called anisotropy matrix obtained in the syn- 
chronous gauge calculations of Ref. [1] . The big difference between our result and that 
of Ref. [1] is the anisotropy matrix of Ref. [1] is independent of time whereas has 
an exponential time dependence. (Clearly, becomes time-independent when the shift 
vector is zero.) Hence, we see explicitly that one can artificially remove or enhance the 
anisotropy coming from , depending on what choice is made for the shift vector. More 
to the point, let us consider the simplified example where only ai, as, and ag are non-zero, 
with ai = as = ag < 0. Then it is easy to see that ^^^JCf decays to zero exponentially with 
time. Hence, as t increases, it will quickly be the case that 



(29) 



The last item for this section is to determine ^°^7y. Since we have seen how to use 
the shift vector to reduce the effect of Qf {x^, t) on ^°)/C^, we will henceforth set both 
and Qf to zero for the rest of this work. Eq. (29) thus becomes the zeroth-order solution 
for the extrinsic curvature. From Eqs. (7), (9), and (29) it follows that 

(0)^. . = (0)^1/3;,^^. ^ h^k^^^ = 51 , (30) 

where hij depends only on a;* and is the so-called "seed" metric mentioned in the intro- 
duction. The remaining field equations only give one more thing, when Eq. (22) is used, 
and that is 

(o)K = -V24^ (o)^-r/4 _ (31) 

To determine the explicit time-dependence of and thus the other zeroth-order quan- 
tities, the lapse function must be specified. For instance if = 1, like in the synchronous 
gauge, then 

(0)^1/3 0, ^4r/3_ Qj.^ condition like K — k{t) must be imposed, which is what 

we do next. 

4. The Long Wavelength Iteration Scheme 

The long wavelength iteration scheme obtains solutions to the Einstein equations by 
assuming spatial gradients of the gravitational and matter fields are small on every spacelike 
slice of the foliation. The solutions contain terms with ever-increasing spatial gradients of 
the "seed" metric obtained in the last section. To have solutions with inhomogeneity on 
shorter and shorter scales, terms with more and more spatial gradients must be included 
[23]. Here we write down solutions with terms containing up to four gradients. 

To be precise, the field variables p, w*, and are expanded as follows (letting Ao{t) = 

(0)^1/3). 

P = Po{t) + p2{t)R + PaR^ + fi4{t)RijR'^ + e^DiD^R + ... , (32) 
Ui = u^{t)DiR + U5{t)RDiR + V5{t)RlDjR + w^{t)R'^DiR]+ 

X5{t)R';DkRi + y^{t)D,DjD^R... , (33) 

and 

Hi = ^'hij + ^'hij + ^'hij + - , (34) 
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where i?^ is the Ricci tensor obtained using the "seed" metric hij, R = h'^^ Rij, Di is the 
covariant derivative compatible with the "seed" hij, 



f2{t)Rij + -{92{t)-f2{t))Rhij 



(35) 
(36) 



and 



(4) 



7y = -4o(t) 



^ (a4(t) - b^{t)) R^Kj + b^{t)RRij + ^ (C4(t) - d^it)) R''^Rkihij + 
d^{t)R]Rki + ^ (e4(0 - h{t) - 94{t)) {DkD'^R) h^j + U{t)DiDjR+ 



g^{t)DkD^Rij\ . 



(37) 



The higher order terms in the expansions above are "small" in the sense that Rij contains 
two spatial gradients. That is, if L is the characteristic scale on which the fields vary, then 
Rij ~ L~'^hij, DiDjRki ~ L~^hijhki, and so on. 

When the K = kit) slicing condition is imposed, the lapse function is no longer freely 
specified, but is obtained from 

N = -^^, (38) 
k{t) V7 ^ ' 

once k{t) is given and the form of 7 is known. The 3+1 field equations also refiect the 
difference, becoming 



ie>TTp = R-K]Kl+k''{t) , 



(39) 

(40) 



1 



-Dib^{\n^--AT,{p-S)5i , (41) 



and 



(42) 



kit) 

Determining from Eq. (38) is unconventional. However, it will be shown to be equivalent 
to the conventional approach (which is to solve Eq. (8)) in the next section. 

The ordinary differential equations that determine the time-dependent coefficients in 
Eqs. (32-34) come about by putting the expansions in Eqs. (32-34) into the 3+1 equations 
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just above, and then rearranging until all "like" terms (i.e., quantities containing hij, Rhij, 
or Rij, etc.) are gathered. One then gets as many equations as coefficients, since at each 
order the tensors Rhij, Rij, etc., are in general linearly independent. We will not list the 
zeroth-order equations and solutions here since they are the same as those presented in 
the last section (except that ^^^K is to be replaced with k{t)). The second-order equations 
for f2{t) and g2{t) are 



■j,2(r-i)/r ffV _ 8 , -2(i+r)/r 



and 



kg'2 = 



rMo(/c) 
2 (2 - 3r) 



where ' = d/dk = k~^d/dt and ^o(^) ^ k~'^/^^ . The fourth-order equations are 

2(3r-2) 2-r 2 k (2 - 3r)^ 



(43) 



(44) 



(45) 



2 Sr _^_2(r+r)/r ^^^^^ ^ ;,,;^/c-^(^+n/r (7^^ + 2^2) - 



r2Ao(A;) 



3r2^o(A;) 



4(2- 3r) , 2(i+2r)/r 
T^Alik) 

1 ' f (1 f'\ ^ 3(2 - r) , 2 , 2(2-3r) 



e^^-^y^ (4 - hf'^) 

3r-2 



16 



rMo(A;) 



^ 3r2/c2^o(A;) 



(/2 - 4^2) + 



4 (3r - 2) 
V^k^Alik) ' 



' 4 (2 - 3r) 2(i^2r)/r 1 ^-2(i+r)/r / <• _ ^ 



and 



where 



V^Ao{k) 



Si^6i + ^%^/^ + i[/2(/2-^2)r • 



(46) 
(47) 
(48) 
(49) 
(50) 

(51) 

(52) 
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The equations for the time dependent coefficients for p and Ji are obtained from Eqs. (70) 
and (71) in Appendix A. Also note that we have still not specified the explicit form for 
kit). 

It is not crucial that we write down here the general solutions to these equations (i.e., 
the sum of the "homogeneous" and "particular" solutions), even though it is straightfor- 
ward to integrate and determine their form [24]. The crucial thing is to notice that Eqs. 
(44), (45), (47), and (49) have only one derivative. In the synchronous gauge approach 
of Ref. [1] the equations all have two derivatives. Therefore, there are "homogeneous" 
solutions appearing in the synchronous gauge calculation that do not appear here. For 
instance, the general solution for g2{k) is 

^2(/c)=7i-^/c^^^-^")/3^ (53) 
whereas the general solution of Ref. [1] is 

To compare these two we must remove the final arbitrariness left in our solutions by 
giving an explicit form for k{t). The appropriate choice is k{t) oc since the zeroth-order 
solution for the extrinsic curvature obtained by Ref. [1] has a trace given by ^^^K oc t~^. 
More generally, since ^o(^) oc then the trace of the extrinsic curvature must decay 

with time if the zeroth-order "expansion factor" Aq is to grow with time. A trace like 
k{t) (X thus results in a power law time dependence for Aq. It now follows that 
the time dependence of the "particular" solutions in Eqs. (53) and (54) is the same 
^^-2(2-3r)/3r^ The synchronous gauge result for g2 has the extra "homogeneous" solution 
72^"^. However, Comer et al. [1] demonstrated a coordinate transformation within the 
synchronous gauge exists whereby this extra term can be removed. The K = k{t) slicing 
condition automatically removes the extra term. It is clear that pure gauge terms will 
automatically be removed at the fourth-order as well. 

Thus, as claimed in the Introduction, there are pure gauge terms appearing in the 
synchronous gauge approach that are not present here. However, we do recover the other 
terms that Ref. [1] obtained. And like the results of Ref. [1] , it can be shown [24] that 
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some of the "homogeneous" solutions (for instance, the 71 term in Eq. (53)), can either 
be absorbed into a redefinition of the "seed" or decay with time for all allowed values 
of r. The other "homogeneous" solutions, as well as the "particular" solutions, decay if 
r < 2/3, that is, when the matter is 'inflationary.' As explained in Ref. [1], this means that 
the metric perturbations will all freeze out once the comoving scale of the Hubble radius 
(which is shrinking like t^~'^/^^) becomes less than the characteristic scale L associated 
with the perturbations. 

When r > 2/3, then the opposite behaviour happens. Instead of decaying with time, 
the metric perturbations will all thaw out, and hence grow, as they enter the Hubble radius. 
In this case, one must keep more and more terms if the solutions are going to be accurate 
[23]. It is for precisely this reason that the K = k{t) condition, and a non-zero shift vector, 
can be of great help. The differential equations for the time-dependent coefficients are 
simpler, and it may be possible to make the expansions converge faster, in the sense that 
a fewer number of iterations would be required. 

5. The Lapse Function Solution 

In this section we will verify that the lapse function used in the previous section does 
indeed satisfy Eq. (8). Recall that we found in Sec. 3 that the spatial gradients of N had 
to be small. With that in mind, N is expanded like 

= No{t) + N2it)R + N^{t)R' + V4it)RiRi + n^{t)DiD'R + ... . (55) 

The equation that determines N is the appropriate form of Eq. (8), setting A'"* = and 
using K = k{t): 

D'DiN - Attn (p + S) - NK{K'j = -k{t) . (56) 

The zeroth-, second-, and fourth-order equations for N are constructed in the same manner 
as the previous section. (For the reader's convenience, the expansions to fourth-order for 
p and KfKj are listed in Appendix A.) But unlike the previous section, the equations for 
No{t), N2{t), etc., are algebraic [24]. 
The solution for N is 



2 k 2-3T . 
T k'' ^ n^Aoik) 



V^k^Ao{ky^ 16A;2 
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kRlR]+ 



2-3r (2-3r 

1/2 - 92) + 



3r2A;4^o(A;) 

2 - 3r 



(/2 - 4^2) + 



2 (2 - 3r) 



kR^+ 



kDiD'R 



(57) 



_6r2fcMo(A;) ' T^k^Al{k)_ 

Note the requirement N > implies k/k'^ > 0. Also notice that Eq. (57) does indeed 
agree with the result obtained from Eq. (38), making use of Eq. (67) from Appendix A as 
well as Eqs. (43-51) for the metric corrections. 

The time-dependence of N can be qualitatively understood by using again k{t) (xt~^ 
and then inserting f2 ~ (72 ~ 

^-2(2-3r)/3r_ jg ggg^ ^Yiat N2/N0 ~ ^-2(2-3r)/3r 
N4^/No r74/iVo ~ n4/No ~ ^-4(2-3r)/3r_ ^hus the corrections to N all decay with time 
when r < 2/3, i.e., the matter is 'inflationary.' Hence, there is no contradiction with 
simultaneously having a monotonic k{t) and a positive lapse function. On the other hand, 
the corrections all grow if F > 2/3. Hence, one must be cautious in the use of the long 
wavelength solutions [23] since the corrections to are not positive-definite. (The second- 
order term is positive at all times when F > 2/3 if the "seed" is such that R < 0.) Certainly, 
N is positive for the amount of time that the second- and fourth-order corrections remain 
smaller than the zeroth-order term. After that, the sign of N will depend on the particular 
value for F and the form of the "seed." 



6. Conclusion 

The main point of the long wavelength iteration scheme is to provide analytical cosmo- 
logical models that contain large scale inhomogeneities. These solutions are perturbative 
and are obtained as expansions in spatial gradients. The fundamental difficulty with this 
scheme is the lack of an absolute notion of space, and thus an absolute notion of spatial gra- 
dient. Any such notion must be introduced by hand. Therefore, it is essential to verify that 
these long wavelength inhomogeneities result from real physics, and not from embedding 
unnecessarily wrinkled spacelike slices into an otherwise homogeneous spacetime. 

We have addressed this issue by applying the long wavelength iteration scheme to 
the 3-1-1 form of the Einstein-Perfect Fluid field equations. We have verified the results 
of Ref. [1] using the K = k{t) slicing condition, and in the process demonstrated that 
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the 3+1 formalism makes it easier to implement the long wavelength iteration scheme. In 
particular, it was seen that a non-zero shift vector can be used to exponentially dampen 

the effect of anisotropy on the extrinsic curvature and also spurious gauge modes that arise 
in the synchronous gauge approach do not show up when the K = k{t) slicing condition 
is invoked. 

Even though we have gone a long way in applying the 3+1 formalism to the long 
wavelength iteration scheme, there are still some questions that need to be considered. For 
instance, we did not solve the 3+1 equations in the canonical way, i.e., as an initial- value 
problem. It is clear that this should be investigated, since ultimately the question of how 
many terms to keep in the expansions should boil down to conditions on the initial data 
set. The problem [25] can be framed around York's procedure for handling the constraints 
[8]. That is, one performs a conformal transformation on the metric, extrinsic curvature, 
and matter variables, and then uses York's covariant decomposition of symmetric tensors 
into their trace-free, vector, and trace parts. There is even a small clue in what has been 
presented here that this is an appropriate way to proceed: the three-metric naturally had a 
conformal factor appear, which was the third-root of the determinant of the three-metric. 
This is precisely the factor suggested by York's procedure. 

Finally, it would also be interesting to see how to construct and use "minimal-strain" 
or "minimal-distortion" shift vectors [10] in the long wavelength iteration scheme. Smarr 
and York [10] show that these vectors are very adept at simplifying the form of the three- 
metric, by separating as completely as possible the purely "kinematical" (in the words of 
Ref. [10]) from the dynamical parts of the three-metric. This is presently being investigated 
[25]. If time-independent "minimal-strain" or "minimal-distortion" shift vectors can be 
constructed, then they can be immediately placed into Eq. (28). 
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Appendix A 

For the convenience of the reader we list here some basic formulas used to build the 
second- and fourth-order equations given in the main text. In all that follows consider that 
the three-metric has the form 

lij = ^""hij + , (58) 

where 

= ^^hij + ^^hij , (59) 

and ^^''7y and ^'^^'fij are the same as Eqs. (36) and (37), respectively, in the main text. 
The first formula is that for the inverse metric 7*-' . It is given by 

where 

X^Aj — {0)^,ik (0)^j7r^, i (0)^,iZ (0) im (0)^,fcn (2)^, (2)^, /^-i \ 

d-f' = ■'-f ^ d-fkl + ^ ^-f ^ ^ ^7 ^ ^■Jmn Ikl ■ (61 j 

The inverse metric is necessary for constructing the connection coefficients, the Ricci tensor 
and scalar, and the extrinsic curvature. 
Letting 

(o)-pi ^ 1 (0) il 
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^'h3i,k + ^'hku 



(0) 



ljk,l 



(62) 



then it can be shown [21] that the connection coefficients are 



(63) 



where 



(64) 



and ^^^Di is the covariant derivative compatible with ^^^^ij {i.e., ^^^Di ^^^^jk = 0). 
The Ricci tensor is obtained from 



(65) 



In terms of the metric corrections given in Eqs. (44) and (45) then 



\g2 (d.D^R - \duD'^R51\ - \f2DkD^ {m - ^RSj ] + 



3/2 ( - i^RfRiSi ] - ^/2i? ( Ri - ^RSi ) + ^ (/2 - 4^2) D.D^R - hRiRt+ 



-{h-g2)RRi 



(66) 



Derivatives of the determinant of the three-metric are obtained from ^7 = ^^'^^S^ij. 
Letting 5 be d/dt, then 



V7 



3Ho + \g2R + ^ (c4 - /2/2) R^Ri + ^ («4 + ^ [/2/2 - ^2^2] ) R^+ 



1 



e^DkD^R 



(67) 



where Hq = AQ{t)/2Ao{t). Notice that it is only g2{t) which enters as a second-order term. 
The extrinsic curvature is 



/2 ( - iRSi) + L - ^/2^2 + I [/2 (/2 - ^2)] ■ ) X 



i? ( i?^' - ^R6i) + (d4 - /2/2) UiRt - iRfRiSi) + h (d.D^R - ^DkD'^ 



g^DuD^ ( i?^ - -R5i 



(68) 



18 



Notice that the correcting pieces are such that we still maintain K = Kl = k(t). It is also 
found that 



1 
3 



{6Ho)^-''' V"'"" 3" ; ■ ^^^^ 

This has corrections beginning with the fourth-order. In the synchronous gauge, there are 
corrections starting with the second-order. However, it is really because K = k{t) holds 
at all orders that our second- and fourth-order equations are simpler than those of the 
synchronous gauge. 

Finally we list the equations — good to fourth-order — for p and Ji: 



167rp = \k\t) + J-i? + ^ [/2 - Ag^] D.D^R - 

6 Aq DAq 



1 

3 



^(/2-^2) + ^^/| 



1 



J_f I fe'(^) f2 

i?2 



(70) 



and 



k{t) ' ' 6 



3^4 - 64 + -^1292 - \ (/2 [/2 - ^2]) 



RDiR+ 



+ -\9^+\ (^4 - /2/2) - -^1292 + ^ (/2 [/2 - ^2]) 



I^iDjR+ 



dA - 72/2 + 45r4 



394 + 



1 
3 



2/4 + 



3 

DkD^DiR 



{d4 - 12/2) 



R'j,DiRj + 



(71) 



From this we see J^iJ^j has six spatial gradients. Hence, UiUj can be ignored in all equations. 
In particular, p ^ p* and ^ (T — 1) pSj. 

Appendix B 



In this appendix we show how to solve Eq. (23) for ^^^JC^. Let 



(72) 



and E°''^ j represent 'a' eigenvectors with eigenvalues A((j), a = 1,2,3, 9; that is. 



(73) 
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with no sum being done over a. Also let [E ■'^) be such that 

j {E-') u ' = and i (^"') al ' = S^Sj . (74) 

a 

Then, the matrix with components Ad'''^ can be diagonlized via 

(^"') w ' = (75) 

Finally, let 

/C« = E'^'/'^^JCi and ^^^X:^' = ^ (S'^) 'K" . (76) 

a 

Therefore, Eq. (23) becomes 

IC" - (L>feA^^) /C" - A(o)/C" = (77) 

and the solutions are thus 

/C« = exp[(i?fciV'= + A(„))t] , (78) 
with only depending on x*. Or, letting 

Gi ^ J]exp(A(.)t) (£;-i) , = , (79) 

we get the solution written in Eq. (28). 
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